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Generalized Grover's searching algorithm for the case in which there 
are multiple marked states is demonstrated on a nuclear magnetic res- 
onance (NMR) quantum computer. The entangled basis states (EPR 
states) are synthesized using the algorithm. 
PACS number(s):03.67 

l.Introdution 

Since the quantum searching algorithm was first proposed by Grover [1], 
several generalizations of the original algorithm have been developed [2]- [4]. 
The generalized algorithm that we will realize can be posed as follows. Let 
A" basis states of a system constitute set D. A function F is denned as 
F : D — ► {0, 1}. The states satisfying F(x) = 1 are defined as marked 
states, which constitute set M with a total of r states. The other states in 
D, satisfying F(x) = 0, constitute set M with a total of A^ — r states. The 
states in M and M have amplitudes fcj and Zj, respectively. A unitary operator 
U transforms a predefined basis state \s > into a superposition denoted as 
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\g(0) >= U\s >. U can be almost any valid quantum mechanical unitary- 
operator. \g(0) > is the initial state for the algorithm. The phase rotation 
of marked states is described by 7 t 7 = J2 X e llF ^ \x >< x\. Obviously, if 
\x >G M, If\x >= e n \x >; if \x >G M, I^\x >= \x >. A composite operator 
is defined as G = —UI^U^. /f is defined as if = I — (1 — e ll3 )\s >< s\, where 
/ denotes unit matrix. The problem is to transform |<?(0) > to a target state 
denoted as \ip t >= J2i^m > by repeating Grover iteration n times. When 
the system lies in \ip t >, measurement yields \ki\ 2 , the probability of the 
system being in marked state \i >. When (5 = n, 7 = n, \s >= |0 >, and U 
is chosen as Walsh-Hadamard (W-H) transform, the generalized algorithm 
becomes the original Grover 's algorithm, where |0 > denotes all qubits in the 
state. 

E.Biham et al have analyzed generalized Grover's algorithm using re- 
cursion equations [3]. Through introducing an ancilla qubit and choosing 
a proper U, Grover proposed a theoretical scheme to synthesize a specified 
quantum superposition on N states in 0(\^N) steps using the algorithm [5]. 
We find that some special superpositions, such as EPR states, can be syn- 
thesized using the algorithm without the ancilla qubit. Generalized Grover's 
algorithm of one marked state has been realized on a two-qubit NMR quan- 
tum computer. G.-L. Long et al realized the algorithm by choosing the phase 
rotation as I t 2 , while the W-H transform is retained [6] . In our previous work, 
we realized the algorithm by replacing the W-H transform by other unitary 
operator, while the n phase rotation was unaltered [7]. In this paper, we 
will realize the generalized algorithm of multiple marked states and synthe- 
size EPR states. The W-H transform is replaced by other unitary operator 
and the phase rotation is chosen as I t 2 or I t 2 . 

2. The generalized Grover's algorithm 

In this section, we will use some results in Ref.[3] to express the principle of 
our experiments. 

n applications of G!} transform \g(0) > into \g(n) >, described by 
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\g(n)>=^k i (n)\i> + , £l i (n)\i>. (1) 
When n — 0, the system lies in the initial state 

10(0) >= U\s > . (2) 

One can find that fcj(0) = Ui S , i G M, and h(0) = Ui S , i G M, where Ui s =< 
i\U\s >. Git transforms the amplitudes kj(n), j G M, to kj(n + 1) =< 
j\GI?\g(n) > and amplitudes lj(n), j G M to lj(n + 1) =< j|G'/ t 7 |5 , (^) > 
The recursion equations describing such iteration are expressed as 



kj{n + 1) = e*(l - e*)^, £ fcj(n)^ + (1 - e l/3 )^ s J2 l i( n ) U ts - e^k,(n), 

i&M i( z M 

(3) 



+ = e*(l-e#)U ja £ Wt/* + (1 - £ /.(n)^* -Z». (4) 

Without loss of generality, we assume that C/ is 7^ 0, « = 1, 2, • • •. Valuables 
k'^n) and I'^n) are defined as 

fcl(n) = (5) 

«(n) = 7^- (6) 

One can easily find that fc^(0) = 1, ^(0) = 1. The weighted averages are 
defined as 

= W E \U is \%(n), (7) 

*>) = WF J2\ a ^\X(n), (8) 

' ieM 

where Wk = Y^ieM \Ui S \ 2 , and Wi = J2 ie M \U is \ 2 . With these variables, the 
recursion equations can be rewritten as 
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k'Jn + 1) = e* 7 (l - e i(3 )W k k'(n) + (1 - e i/J )W,Z» - e^k'An), (9) 



Z^(n + 1) = e i7 (l - e if> )W k k'(n) + (1 - e i/3 ) W» - I'^n). (10) 

By averaging over all the marked states in Eq. (9) and over all the unmarked 
states in Eq.(10), we find the two recursion equations for k'(n) and I' in) can 
be expressed as 



k'(n + 1) = e* 7 (l - e^)W k k'(n) + (1 - e i/3 ) W(n) - e l ^k'(n), (11) 

7{n + 1) = e* 7 (l - e i/3 )W fe F(n) + (1 - e i/J )WjZ>) - F (n). (12) 
Subtracting Eq.(ll) from Eq.(9), and Eq.(12) from Eq.(10), one finds that 

k'jin + 1) - F(n + 1) = -e l7 (^.(n) - F(n)), (13) 

/;.(n + 1) - 7(n + 1) = -(/;.(n) - 7 (n)). (14) 

Noting that fcj(0) = 1, Z|(0) = 1, we find F(0) = 1, f (0) = 1. Using Eqs.(13) 
and (14), we obtained 

k' l (n)=¥(n), (15) 

l'i(n) =7(n), (16) 

where the subscript j in Eqs.(13) and (14) has been replaced by %. 

From the discussion above, for any U, k'(n) and l'{n) can be solved from 
Eqs. (11) and (12). Using Eqs.(15),(16),(5),and (6), we can obtain the 
explicit expressions for fcj(n) and li(n). Eqs. (11) and (12) can be rewritten 

as 

U(n+1) j' A {7(n) J' (17) 
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where 



A= ( e*(l - e V)W k - e* (1 - e^)W, \ 

Eq.(17) shows that fc'(n) and /'(n) are dependent on 7, /3, and If U is 

replaced by a different U' where \U' is \ = \Ui s \, the analysic forms of k'(n) and 
l'(n) are unaltered. If kin) approaches when n = n , the system lies in 

\1>t>= E^K)K>. (19) 

Generally, the state is not the equally weighted superposition of marked 
states. It is related to Ui S . By choosing proper U, some superpositions can 
be synthesized using generalized Grover's algorithm. We will solve Eq.(17) 
in the following section. 



3. Experimental scheme 

Our experiments use a sample of Carbon- 13 labelled chloroform dissolved in 
d6-acetone. Data are taken at room temperature with a Bruker DRX 500 
MHz spectrometer. The resonance frequencies v\ = 125.76 MHz for 13 C, and 
v 2 = 500.13 MHz for 1 H. The coupling constant J is measured to be 215 
Hz. If the magnetic field is along i-axis, and let h — 1, the Hamitonian of 
this system is described by 

H = -2-kvJI - 2nu 2 I 2 z + 2irJlll*, (20) 

where I^(k = 1,2) are the matrices for ^-component of the angular momen- 
tum of the spins [8]. In the rotating frame of spin k, the evolution caused 
by a radio-frequency (rf) pulse on resonance along x or — y-axis is denoted as 
Xk(<Pk) — e t(pkI * or Yj~{— (pk) = e~ tlpkI y , where ipk = B\^ k t p with k specifying 
the affected spin. B±, 7^ and t p represent the strength of rf pulse, gyromag- 
netic ratio and the width of rf pulse, respectively. The pulse used above is 
denoted as [ip]* or [— <p]y. The coupled-spin evolution is denoted as 

[t] = e- <2wJJ « 1J *, (21) 
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where t is evolution time. The predefined pseudo-pure state 



\s >= | T>i I T>s 



o 



Vo ) 



(22) 



is prepared by using spatial averaging [9], where | \> k denotes the state of 
spin k. For convenience, the notation | |>x | |> 2 is simplified as | Tt>- 
The basis states are arrayed as | Tt>> I tl>> I IT> 5 I ll>- U is chosen as 
U = Yi((pi)Y 2 ((p2) represented as 



U = 



( C\C2 C\S2 S\C2 S\S2 \ 

-CiS 2 CiC 2 SiS 2 SiC 2 

-SiC 2 -SiS 2 CiC 2 CiS 2 

V S1S2 -sic 2 -cis 2 cic 2 / 



(23) 



where = cos({p k /2), s k = sin((p k /2). When | |f> and | H> are the two 
marked states, U can be chosen as U = Y 1 (^)Y 2 (^) described by 



/ 1 



U 



One can find that W k 
represented as 



1 

~ 2 
Wi 



1 1 
-1 1 -1 
-1 -1 1 



1 \ 
1 
1 



(24) 



V 1 -1 -1 1 J 



When 7 = = -f , 77 and If can be 



[ 14 



/ -i 





\ 

1 
1 



V 



(25) 



/ -i \ 

10 

10 

V 001/ 



(26) 
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Using the values of f3, 7, W k , and W h Eq.(18) is expressed by 

A= 1 -( 1+i (27) 
2\l-i i-1 J K J 

In order to obtain explicit expressions for k'(n) and l'(n), we introduce the 
diagonal matrix represented as 

A D = S-'AS = ( A + A ° ) (28) 

The eigenvalues of matrix A are the solutions of det(A — XI) = 0. They are 
expressed as A + = e l ? , A_ = e l ~%~ . S and S^ 1 are expressed as 



(29) 



1 


1 


V3-1 


v^+1 


1+i 


1+i 


1 v^+1 


1+i 








1+i 


V 2V3 


2V3 



The solution of Eq.(17) can be expressed as 

(KH(S)< 

where A n = SA^S' 1 , expressed as 



A n J_ ( (VS + l)e^ + (VS - l)e^ (l + i)(e^-e^) \ 

2y/3\ (l-i)(e^- e ^) (y3-l)e^ + (V3 + l)e^ J ' 

(32) 

When n = 1, we obtain that k'(l) = v^e 1 *, and /'(l) = 0. Using C/n = |, 
C/41 = |, we obtain that fci(l) = e % ^/y/2, and £4(1) = e l ^/y/2. The system 
lies in state 

l^i >= (I TT> +1 !!>)e^/v / 2. (33) 
The overall phase can be ignored. 
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If U is chosen as U = Yi(-§ )Y 2 (^), Eq.(32) is unaltered. We also obtain 
fc'(l) = y/2e tJ *, and /'(l) = 0. Noting that Un = \, and U 41 = — |, we obtain 
that fci(l) = e l %/y/2, and fc 4 (l) = —e l ^/V2. The system lies state 



life >= (I TT> -I U>KVv2. (34) 

Considering the experimental convenience, if the marked states are | ||> 
and | |t>) we choose I? as I] 3 = «/ 14 2 , where 7 = f • In matrix notation, I] 3 
is represented as 



T 2 



/ 1 



Vo 



\ 




1 / 



(35) 



(3 is changed to | to satisfy phase matching [10]. Is is described by 



T 2 — 



I % \ 
10 
10 

V o o o i j 

When 7 = /3 = § , and [/ = Yi(|)y 2 (|), or [/ = 
of Eq.(17) can be obtained by replacing % in Eqs.(27)-(32) by —%. We obtain 
F(l) = V2e-*5, and = 0. When ?7 = Yi(f )F 2 (f ), we obtain fc 2 (l) = 
— e~ l ^/y/2, and fc 3 (l) = — e~ l ^/V2, using C/ 2 i = — \, and C/ 3 i 
system lies in state 



(36) 



y"i(-f )l2(f), the solution 



-\. The 



life >= "(I U> +1 iT>)e-'/V2. 



(37) 



Similarly, when U = F^-f^d ), we obtain fc 2 (l) = -e^/v 7 ^, and M 1 ) 



2 7 v 2 y 

e _l ?/-\/2, using C/ 2 i = — |, and C/ 3 i = |. The system lies in state 



life >= "(I U> -I IT>)e-^/V2. (38) 

| -01 >, |^2 >) 1^3 > an d |ife > are the four EPR states. They are very 
useful in quantum information and have been implemented in experiments 



8 



[11] [12]. Based on the discussion above, they can be synthesized by general- 
ized Grover's algorithm. Other entangled states can be obtained by choosing 
other U . For example, if U is chosen as 



/lli i \ 
-1 1 -i i 
till 
\ -i i -1 1 / 



(39) 



and 7 = (5 — — |, entangled state (| |T> ~A ||>)e^/v^2 is obtained after 
one iteration. The target states, such as \ip\ > and \ip 2 >, can also be 
obtained by matrix multiplication. If replacing n in Eq.(32) by n + 3, one 
finds A n+3 = iA n . This fact shows that k'(n) and l'(n) both have a period 
of 3. 



4. Experimental procedure 

The equilibrium density matrix can be represented as 



Peg 



7T 



(40) 

1,2 _ 



The rf and gradient pulse sequence [a\1 — [grad] z — [tt/A]]. — 1/4J 
1/4 J— [— 7r]J.' 2 — [— 7r/4]^ — [grad] z transforms the system from the equilibrium 
state into the state represented as 



p s = 11/2 + PJ2 + III 



( 3 



Vo 



000 
-100 
0-10 
00-1 



(41) 



which can be used as the pseudo-pure state | ft> [13]. a = arccos (71/272), 
[grad] z denotes gradient pulse along z- axis, and the symbol 1/4 J means 
that the system evolutes under H described as Eq.(20) for 1/4J time when 
pulses are closed. The pulses are applied from left to right. [^}l' 2 de- 
notes a nonselective pulse (hard pulse). The evolution caused by the pulse 
sequence 1/4 J — [vr]^ 2 — 1/4 J — [— n]l' 2 is equivalent to the coupled-spin 
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evolution [1/2 J] described in Eq.(21) [14]. n pulses are applied in pairs 
each of which take opposite phases in order to reduce the error accumu- 
lation caused by imperfect calibration of [tt] pulses [15]. U is realized by 

7T 

[±7r/2]i — [7r/2] 2 , corresponding to Yi(±j)Y 2 (^), respectively. Because 1 2 \ — 
i/ 14 2 , I 2 \ and J 14 2 can be realized by the same sequence 1/4 J — [7r]*' 2 — 
1/4J — [— 7i]l' 2 . By modifying the pulses used in Refs.[6][16], we realize I s 2 
by 1/8 J - K> 2 - 1/8 J - [-vr]^ - [-tt/2]^ - - [vr/2]^ 2 , and i} 

by 15/8 J - K< 2 - 15/8 J - [-tt] ^ - [-tt/2]}' 2 - [vr/4]^ 2 - [vr/2]^ 2 . When 
U = Yi(|)y 2 (|), (-Uh^U^I^)U transforms the pseudo-pure state | tt> 

2L 2L 

into state [^i >, and (—UIsWl 23 )U transforms | ||> into >, where 
U transforms | j"f> into the initial state \g{0) >, and () indicates Grover 
iteration. When U = Fi(-f )F 2 (f ), (-Ul7*U*I^)U transforms | |T> into 
I -02 >, and {—UIsU^I 2 \)U transforms | ft> into \ip4 >. The results are ex- 
pressed by density matrixes. For example, the density matrix corresponding 
to |^i > is represented as 



P1 = mil - im + iiit) = 



( 0.25 0.5 \ 

-0.25 

-0.25 

V 0.5 0.25 j 



(42) 



which is equivalent to >< A readout pulse [vr/2] 2 transforms p\ into 
pi r represented as 



/ 



Plr 



V 



-1 

1 -1 
1 -1 



1 1 \ 

-1 -1 

1 

1 / 



(43) 



The information on matrix elements (1,3) and (2,4) in Eq.(43) can be directly 
obtained in the carbon spectrum, and the information on elements (1,2) and 
(3,4) can be directly obtained in the proton spectrum. Similarly, when the 
system lies in \ip 2 >, \ip 3 >, or >, the readout pulse [tt/2] 2 transforms the 
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system into the state represented as 




1 



/ — 1 —1 —1 \ 
-10 11 
-110 1 



(44) 



V -l i i o j 



i 



( 
i 



i i -l \ 

1 -1 

1 -1 
-1-10/ 



(45) 




1 



V -i 



or 



4 



1 



/ 1 -1 1 \ 
10-11 
-1 -1 -1 



(46) 



Vi i -l o J 



Through observing the matrix elements (1,3), (2,4), (1,2) and (3,4) in Eqs.(43)- 
(46), one can distinguish the four EPR states. 



In experiments, for each target state, the carbon spectrum and proton spec- 
trum are recorded in two experiments. For different target states, carbon 
spectra or proton spectra are recorded in an identical fashion. Because the 
absolute phase of an NMR signal is not meaningful, we must use reference sig- 
nals to adjust carbon spectra and proton spectra so that the relative phases 
of the signals are meaningful [17]. When the system lies in the pseudo-pure 
state described as Eq.(41), the readout pulses [ir/2] y and [n/2] y transform it 
into states represented as 



5. Results 



Psrl — ^ 



1 



/ 1 -2 \ 
0-100 
-2010 

V o o o -l j 



(47) 



11 



and 




1 



/ 1 -2 \ 
-2100 
0-10 
V -1 j 



(48) 



respectively. In the carbon spectrum or proton spectrum, there is only one 
MNR peak corresponding to element (1,3) in p srl or to element (1,2) in p sr2 . 
Through calibrating the phases of the two signals, the two peaks are adjusted 
into absorbtion shapes which are shown as Fig. 1(a) a for carbon spectrum and 
Fig. 1(b) for proton spectrum. The two signals are used as reference signals of 
which phases are recorded to calibrate the phases of signals in other carbon 
spectra and proton spectra, respectively. One should note that the minus 
elements in Eq.(47) and Eq.(48) are corresponding to the positive peaks in 
Fig. 1(a) and Fig. 1(b). 

We implement generalized Grover's algorithm starting with the initial 
state \g(0) >. Git transforms \g(0) > into one of EPR states. If no readout 
pulse is applied, the amplitudes of peaks is so small that they can be ignored. 
By applying the spin-selective readout pulse [vr/2]^, we obtain the carbon 
spectra as shown in Figs. 2 (a), (b), (c), and (d), and the proton spectra as 
shown in Figs. 3 (a), (b), (c), and (d). Fig. 2 (a) and Fig. 3 (a) are corresponding 
to >, Fig.2(b) and Fig.3(b) to \ip 2 >, Fig.2(c) and Fig.3(c) to \i/j 3 >, 
and Fig. 2(d) and Fig. 3(d) to >. In Fig. 2(a), for example, the right 
and left peaks are corresponding to the matrix elements (1,3) and (2,4) in 
Eq.(43), respectively. Similarly, in Fig. 3 (a), the two peaks are corresponding 
to the matrix elements (1,2) and (3,4) in Eq.(43). The phases of the signals 
corroborate the synthesis of EPR states. 

6. Conclusion 

In experiments, we demonstrate generalized Grover's algorithm of 2 marked 
states. The results show that generalized Grover's algorithm is efficient for 
the case of N/2 marked states. The original Grover's algorithm, however, 
does not work for this case. EPR states are synthesized using the algorithm. 
For the case of multiple marked states, the final signal is an average over all 
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the marked states. It is difficult or impossible to deduce anything about indi- 
vidual marked states from the ensemble average [18]. In our work, however, 
the generalized Grover's algorithm is viewed as a technique for synthesizing 
a particular kind of superposition of marked states [5]. The superposition is 
described as a density matrix of which elements can be obtained in the NMR 
spectrum by readout pulses [19]. 
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Figure Captions 



1. The carbon spectrum (Fig. 1(a)) obtained through selective readout 
pulse for 13 C [7r/2]* and the proton spectrum (Fig.2(b)) obtained through 
selective readout pulse for 1 H [n/2]y when the two-spin system lies in 
pseudo-pure state | Tt>- The two peaks are adjusted into absorbtion 
shapes. The two signals are used as reference signals to adjust other 
spectra. 

2. Carbon spectra obtained through [n/2]y after EPR states are synthe- 
sized. Figs. 2 (a), (b), (c) and (d) are corresponding to states (| ft> 
+ | ||>)/v / 2, (| TT> -I ||>)/v / 2, (I T!> +| IT>)/v / 2, and (| T!> 

— | ||>)/v / 2, respectively. 

3. Proton spectra obtained through [vr/2]^ after EPR states are synthe- 
sized. Figs. 3 (a), (b), (c), and (d) are corresponding to states (| |f> 
+ | ||>)/ v / 2, (| TT> -I ||>)/v / 2, (I T!> +| IT>)/v / 2, and (| U> 

— | ||>)/v / 2, respectively. 

[Figure 1 about here.] 
[Figure 2 about here.] 
[Figure 3 about here.] 
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